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Abstract. In a previous work [18l [19] it is investigated, by means of computational 
QJ simulations, shapes of nonconvex bodies that maximize resistance to its motion on a 

rarefied médium, considering that bodies are moving forward and at the same time slowly 
rotating. Here the previous results are improved: we obtain a two-dimensional geometric 
C\J shape that confers to the body a resistance very close to the supremum value (i? = 

1.4965 < 1.5). 

CL| Resumo 

i^H Um corpo bidimensional, apresentando um ligeiro movimento rotacional, desloca- 

is se num meio rarefeito de partículas que colidem com ele de uma forma perfeitamente 

^ elástica. Em investigações que os dois primeiros autores realizaram anteriormente |181 119j , 

I— I procuraram-se formas de corpos que maximizassem a força de travagem do meio ao seu 

movimento. Dando continuidade a esse estudo, encetam-se agora novas investigações que 
^ culminam num resultado que representa um grande avanço qualitativo relativamente aos 

ly-j então alcançados. Esse resultado, que agora se apresenta, consiste numa forma bidimen- 

I sional que confere ao corpo uma resistência muito próxima do seu limite teórico. Mas o 

IO seu interesse não se fica pela maximização da resistência newtoniana; atendendo às suas 

características, apontam-se ainda outros domínios de aplicação onde se pensa poder vir 
a revelar-se de grande utilidade. Tendo a forma óptima encontrada resultado de estudos 
numéricos, é objecto de um estudo adicional de natureza analítica, onde se demonstram 
algumas propriedades importantes que explicam em grande parte o seu virtuosismo. 

> 
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1 Introdução 

Uma área de investigação da Matemática contemporânea ocupa-se com a procura de formas 
de corpos, dentro de classes predefinidas, que permitam minimizar ou maximizar a resistência 
a que ficam sujeitos quando se desloquem em meios rarefeitos. O primeiro problema desta 
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natureza remonta já à década de 80 do século XVII, altura em que Isaac Newton estudou, em 
|14) . um problema de resistência mínima para uma classe específica de corpos convexos, que 
se deslocassem em meios de partículas infinitesimais, de tal modo rarefeitos que fosse possível 
negligenciar qualquer interacção entre as partículas, e que a interacção destas com o corpo 
pudesse ser descrita por colisões perfeitamente elásticas. Mais recentemente temos assistido 
a desenvolvimentos importantes nesta área com a generalização do estudo a novas classes de 
corpos e a meios com características menos restritivas: problemas de resistência em corpos 
não simétricos [31 Hl |5l [TTl [12] , em corpos não convexos de colisões singulares [21 |5l |6l [13] e de 
colisões múltiplas [13] HHl HZ] > corpos de superfície desenrolável [13] , considerando colisões com 
atrito [9] e em meios de temperatura positiva [20] . Mas, quase invariavelmente, os resultados 
que têm vindo a ser publicados têm dado especial atenção a classes de corpos convexos. 

A convexidade de um corpo é uma condição suficiente para que a resistência seja função 
unicamente de colisões singulares — todas as partículas colidem uma só vez com o corpo. Esse 
atributo permite reduzir consideravelmente a complexidade dos problemas tratados. Mesmo 
os vários estudos sobre classes de corpos não convexos que têm surgido, especialmente na 
última década, assentam quase sempre em condições que garantem um único impacto por 
partícula — [21 [3 [H [13] . Só muito recentemente, começaram a surgir alguns estudos prevendo 
múltiplas reflexões, como é o caso dos trabalhos de Plakhov |15| [TUI IT7j . 

Na classe de corpos convexos, o problema reduz- se normalmente à minimização da fun- 
cional de Newton — uma fórmula analítica para o valor da resistência, proposta por Isaac 
Newton [2]. Mas no contexto de corpos não convexos não existe qualquer fórmula simples 
conhecida para o cálculo da resistência. Ainda que seja extremamente complexo, em geral, 
tratar analiticamente problemas de múltiplas colisões, para alguns problemas de minimização 
específicos a tarefa não se tem revelado particularmente difícil, existindo inclusive já alguns 
resultados disponíveis |151I16] . Se, pelo contrário, considerarmos o problema de maximização, 
então a solução chega mesmo a ser trivial — para qualquer dimensão, basta que a parte frontal 
do corpo seja ortogonal à direcção do movimento. 

E se o corpo exibir, para além do seu movimento de translação, um ligeiro movimento 
rotacional? Quando pensamos neste tipo de problemas temos, por exemplo, em mente satélites 
artificiais, de órbitas relativamente baixas, que não disponham de qualquer sistema de controlo 
que estabilize a sua orientação, ou outros engenhos em condições semelhantes. Nessa situação, 
imaginamos que, ao longo do seu percurso, o engenho rode lentamente sobre si próprio. 

O problema de minimização da resistência média em corpos rotativos não convexos foi já 
estudado para o caso bidimensional [T7]: demonstrou-se que a redução da resistência que é 
possível obter em relação ao caso convexo não ultrapassa os 1.22%. Por sua vez, o problema 
de maximização da resistência média de corpos em rotação está longe de ser trivial, contrari- 
amente ao que se passa quando se trata de movimento puramente translacional. Foi, por isso, 
esta classe de problemas o objecto de estudo do trabalho realizado pelo autores em [TH] [l9] : 
investigaram-se formas de corpos não convexos que maximizassem a resistência que os mes- 
mos teriam que enfrentar quando se deslocassem em meios rarefeitos e, simultaneamente, 
exibissem um ligeiro movimento rotacional. Com o estudo numérico que se levou a cabo, 
foram encontradas várias formas geométricas que conferiam aos corpos valores de resistência 
bastante interessantes; mas foi em investigações posteriores, desenvolvidas na continuação 
desse trabalho, que os autores conseguiram chegar ao melhor dos resultados: uma forma bidi- 
mensional que confere ao corpo uma resistência muito próxima do seu limite máximo teórico. 
E este último resultado que agora aqui se apresenta. 

A apresentação do trabalho encontra-se organizada como se segue. Em ^começamos por 
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definir, para o caso bidimensional, o problema de maximização, objecto do presente estudo. 
Depois, em ^ descrevemos o estudo numérico que foi realizado no encalço do corpo de re- 
sistência máxima e apresentamos o principal resultado original deste trabalho: uma forma 
bidimensional que maximiza a resistência newtoniana. A forma bidimensional encontrada é 
então objecto de um estudo em í|4] onde se demonstram algumas propriedades importantes 
que ajudam a explicar o valor de resistência que apresenta. Em ^ apresentamos um estudo 
exploratório sobre outras possíveis aplicações do nosso resultado. Por fim, em f|6] apresenta- 
mos as principais conclusões do presente trabalho e incluímos alguns apontamentos sobre as 
possíveis direcções do trabalho a realizar futuramente. 



2 Definição do problema para o caso bidimensional 

Considere-se um disco em rotação lenta e uniforme, deslocando-se numa direcção paralela ao 
seu plano. Denotemos o disco de râio T por e a sua fronteira por dCr- Retiremos então 
pequenas porções do disco ao longo de todo o seu perímetro, numa vizinhança e de dCr, 
com e G M_|_ de valor arbitrariamente pequeno quando comparado com o valor de r. Ficamos 
assim com um novo corpo B definido por um subconjunto de Cr e caracterizado por uma certa 
rugosidade ao longo de todo o seu perímetro. A questão essencial que se coloca é a seguinte: 
até quanto pode ser aumentada a resistência de um corpo B7 Mais do que conhecermos o 
valor absoluto dessa resistência, estamos sobretudo interessados em saber qual o ganho que se 
consegue obter em relação à resistência do corpo liso (contorno perfeitamente circular, neste 
caso), ou seja, conhecer o valor normalizado 

^ Resistência(i?) 

Resistência(Cr-) 

É possível, desde logo, conhecermos alguns valores de referência importantes para a resistência 
normalizada: R{Cr) = 1 e o valor da resistência R{B) terá que se situar entre 0.9878 (|il7j) 
e 1.5. O valor 1.5 será hipoteticamente atingido se todas as partículas forem reflectidas pelo 
corpo com uma velocidade v"*" (velocidade com que as partículas se afastam definitivamente do 
corpo) oposta à velocidade de incidência v (velocidade com que as partícula atingem o corpo 
pela primeira vez), = — v, situação em que é transmitida ao corpo a máxima quantidade de 
movimento. E ainda possível conhecermos o valor da resistência de alguns corpos elementares 
do tipo B. E o caso, por exemplo, de discos com o contorno integralmente formado por 
reentrâncias rectangulares arbitrariamente pequenas ou com a forma de triângulos rectângulos 
isósceles. Tal como demonstrado em jlHl HU], esses corpos têm associado uma resistência, 
respectivamente, de i? = 1.25 e R = V2. 

Para além de ser definido numa vizinhança e interior da fronteira do disco Cr, assume-se 
que o corpo que se pretende maximizar é um conjunto i3 G conexo e limitado, e com 
fronteira dB seccionalmente suave. Considere-se então um bilhar em M? \ B. Uma partícula 
infinitesimal move-se livremente, até que, ao colidir com o corpo B, sofre várias reflexões (uma 
no mínimo) em pontos regulares da sua fronteira dB, acabando por retomar um movimento 
livre que a afasta definitivamente do corpo. Represente-se por convi? o conjunto definido pelo 
invólucro convexo de B. A partícula intercepta por duas vezes o contorno 9(convi3): quando 
entra no conjunto convB e no momento em que o abandona. Considere-se L = \d{convB)\ o 
comprimento total da curva d{convB) e represente-se por v e v"*" a velocidade da partícula no 
primeiro e segundo momentos de intercepção, e x e x"*" os respectivos pontos onde ocorrem. 
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Denote-se ainda por ip e ip~^ os ângulos que os vectores — v e v"*" fazem com o vector normal à 
secção da curva d{comíB) (direccionado para fora) entre os pontos x e . Serão positivos se 
forem definidos no sentido anti-horário a partir do vector normal, e negativos caso contrário. 
Com estas definições, tanto como ip'^ tomam valores no intervalo [— 7r/2, 7r/2]. 

Representando as cavidades que caracterizam o contorno de B por subconjuntos 
r2i,r22,..., que no seu todo formam o conjunto comiB \B, a, resistência normalizada do 
corpo B (equação ([TJ ) assume a seguinte forma (cf . [ISl EH] ) : 

sendo Lq = |9(convi?) n dB\ o comprimento da parte convexa do contorno dB, Li = 
|9(convi?) n com i = 1,2,..., o tamanho da abertura da cavidade Oj, e -R(ííj) a re- 
sistência normalizada da cavidade ííj, em relação a um segmento liso de tamanho unitário, 
com 

i?($7j) = - / / (^1 + cos [ip^{x,(p) — ip)) cos ífdipdx. (3) 

o J-1/2 J-tt/2 

A função deve ser vista como o ângulo de saída duma partícula que interage com uma 
cavidade (li de abertura de tamanho unitário, centrada no eixo dos x (ver ilustração da 
Figura[T| , e com a forma de — relacionando-se com esta por um factor de proporcionalidade 
l/U. 




Figura 1: Exemplo de trajectória que uma partícula descreve no interior de uma cavidade 
normalizada. 



Da equação (|2]) percebemos que a resistência de um corpo B pode ser vista como uma 
média ponderada (^^ Li/ L = 1) das resistências das cavidades individuais que caracterizam 
toda a sua fronteira (incluindo "cavidades" lisas), multiplicada por um factor que relaciona os 
perímetros dos corpos convB e Cr- Assim, maximizarmos a resistência do corpo B equivale 
a maximizarmos o perímetro de convi? (\d{convB)\ < \dCr\) e as resistências individuais das 
cavidades íí,. 

Encontrada a forma óptima íí*, que maximize a funcional ([s]), o corpo B de resistência 
máxima será aquele cuja fronteira seja formada unicamente pela concatenação de pequenas 
cavidades com essa forma. Podemos então restringir o nosso problema à subclasse de corpos 
B que tenham a sua fronteira integralmente preenchida por cavidades iguais, e com isso 
admitir, sem qualquer perda de generalidade, que cada cavidade ííj ocupa o lugar de um arco 
de círculo de tamanho e <^ r. Como Lj = 2rsin(e/2r), a razão entre os perímetros assume o 
valor 

|9(conv5)| sin(e/2r) ^ i^/r)'^ 
\dC^\ e/2r ^ 2^' 
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ou seja, dado um corpo B de fronteira formada por cavidades com a forma Í7, de ([2]) e (|4j) 
concluímos que a resistência total do corpo será igual à resistência da cavidade individual Í7, 
menos uma pequena fracção desse valor, que será negligenciável quando e <C r, 

R{B) « R{n) - (5) 

Assim, as nossas pesquisas têm como objectivo encontrar formas de cavidades íl que 
maximizem o valor da funcional ([s]), cujo supremo sabemos situar-se no intervalo 

1 < supí^i?(íí) < 1.5, (6) 

como facilmente se comprova usando ([3|: se for um segmento liso, ip'^{x,ip) = —(p e 
= ^ f^y^ f^^J^^:^ {1 + cos {2íp)) cos ip díf dx = 1; nas condições de resistência máxima 

ip~^{x,íp) = ip, logo < J^/^i^^ J^^^i^2 cos <p dip dx = 1.5. 



3 Estudo numérico do problema 

Na classe de problemas que estamos a considerar, apenas para algumas formas muito 
elementares é possível desvendar a fórmula analítica da sua resistência (|3]), como são disso 
exemplo as formas rectangulares e triangulares anteriormente referenciadas. Para formas 
um pouco mais elaboradas, o cálculo analítico torna-se rapidamente demasiado complexo, 
senão impossível, dada a grande dificuldade em conhecermos a função 99"*" : [—1/2,1/2] x 
[— 7r/2,7r/2] [— 7r/2, 7r/2], que, como sabemos, está intimamente relacionada com o formato 
da cavidade O. Assim, o recurso à computação numérica surge como a abordagem natural e 
inevitável para se poder investigar essa classe de problemas. 

Desenvolveram-se modelos computacionais que simulam a dinâmica de bilhar no interior 
de cada uma das formas estudadas. Os algoritmos de construção desses modelos, bem como 
os responsáveis pelo cálculo numérico da resistência associada, foram implementados usando 
a linguagem de programação C, dado o esforço computacional envolvido (a linguagem C foi 
criada em 1972 por Dennis Ritchie; para o seu estudo sugerimos, entre a extensa documentação 
disponível, aquele que é o livro de referência da linguagem, da autoria de Brian Kernighan e do 
próprio Dennis Ritchie, |lQj). A eficiência do código objecto, gerado pelos compiladores de C, 
permitiu que a aproximação numérica de ([s]) fosse realizada com um número suficientemente 
elevado de subdivisões dos intervalos de integração — entre algumas centenas e vários milhares 
(até 5000). Os resultados foram, por isso, obtidos com uma precisão que atingiu em alguns 
dos casos 10^^. Essa precisão foi controlada por observação da diferença entre aproximações 
sucessivas da resistência R que se iam obtendo com o aumento do número de subdivisões. 

Para a maximização da resistência dos modelos idealizados, usaram-se os algoritmos de 
optimização global da toolhox ^^Genetic Algorithm and Direct Search" (versão 2.0.1 (R2006a), 
documentada em [1]), uma colecção de funções que estende as capacidades de optimização do 
sistema de computação numérica MATLAB. A opção pelos métodos Genéticos e de procura 
Directa deveu-se essencialmente ao facto de os mesmos não requererem qualquer informação 
acerca do gradiente da função objectivo nem de derivadas de ordem superior — como a forma 
analítica da função resistência é em geral desconhecida (dado que depende de ip~^{x, ip)), esse 
tipo de informação, caso fosse necessária, teria que ser obtida por aproximação numérica, 
algo que dificultaria imenso o processo de optimização. O sistema de computação MATLAB 
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(versão 7.2 (R2006a)) foi também escolhido por dispor de funcionalidades que permitiram que 
se usasse para função objectivo a subrotina compilada em C de cálculo da resistência, bem 
como a função 92^(3;, ip) por si invocada. 

3.1 "Dupla Parábola": uma forma bidimensional que maximiza a re- 
sistência 

No estudo numérico que os autores realizaram em |18| [TÕ] procuraram-se formas íí/ definidas 
por funções / : [—1/2, 1/2] contínuas e seccionalmente diferenciáveis: 

íí/ = {(x, y) : -1/2 <x< 1/2, O < y < f{x)} . (7) 

Iniciou-se a procura da resistência máxima na classe de funções contínuas / com derivada /' 
seccionalmente constante, alargando-se depois o estudo a classes de funções com a segunda 
derivada /" seccionalmente constante. No primeiro dos casos o contorno de ííj é uma linha 
poligonal, e no segundo, uma curva composta por arcos de parábolas. Não se tendo con- 
seguido com as formas íí/ superar o valor de resistência R = 1.44772, decidimos, neste novo 
estudo, estender a procura a formas diferentes das consideradas em ([7|. Estudámos formas 
definidas por funções de y em x do seguinte modo: 

09 = {(x, y) : O < y < /i, -g{y) <x< g{y)} , (8) 

onde h > O e g : [O, h] ^ M.^ é uma função contínua com g{0) = 1/2 e g{h) = 0. 

O novo problema de resistência máxima por nós estudado pode então ser formulado da 
seguinte forma: 

Encontrar supgi?(í79) nas funções g : [0,h] — > contínuas e seccionalmente 
diferenciáveis, tais que g{0) = 1/2 e g{h) = O, com h > 0. 

À semelhança do estudo que se fez para os conjuntos ííj, na procura das formas íí^ 
consideraram-se funções g seccionalmente lineares e seccionalmente quadráticas. Se na classes 
das funções lineares não se conseguiu qualquer ganho de resistência relativamente aos resul- 
tados obtidos para os conjuntos íí/, já nas funções quadráticas os resultados superaram as 
melhores expectativas: encontrou-se uma forma de cavidade íí^ que apresenta uma resistência 
R = 1.4965, um valor já muito próximo do seu limite teórico de 1.5. Este é seguramente um 
resultado muito interessante. Efectuaram-se ainda alguns testes com funções polinomiais de 
ordem superior ou descritas por secções cónicas específicas mas, não se tendo verificado qual- 
quer ganho adicional na maximização da resistência, optou-se por não se reportar os respec- 
tivos resultados. Segue-se então a descrição do melhor resultado que se obteve, encontrado 
na classe das funções x = ±g{y) quadráticas. 

Estudou-se o valor da resistência de conjuntos íí^, tal como definidos em (|8|, na classe de 
funções quadráticas, 

9hAy) = + Py + 1/2, para O <y <h, 

onde h > O e a = ^^^^2^^"^ (dado que gh,/3{h) = 0). Na optimização da curva, fizeram-se 
variar os dois parâmetros de configuração da função: h, a altura da curva dQ^, e /?, o seu 
declive na origem (g'(0)). Nesta classe de funções os algoritmos de optimização convergiram 
rapidamente para um resultado muito interessante: a resistência máxima foi atingida com 
h = 1.4142 e (3 = 0.0000, e assumiu o valor R = 1.4965, ou seja, um valor 49.65% acima da 
resistência da forma rectilínea. E efectivamente um resultado muito importante: 
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i. representa um ganho considerável no valor da resistência, relativamente ao melhor re- 
sultado anterior (obtido em [IHl[l9]), que se situava 44.77% acima do valor de referência; 



ii. o conjunto correspondente tem uma forma bastante mais simples do que a do con- 
junto Qf associada ao melhor resultado anterior, uma vez que é formada por dois arcos 
de parábolas simétricos, quando a anterior era constituída por catorze desses arcos; 

iii. esse novo valor de resistência está já muito próximo do seu limite máximo teórico que, 
como se sabe, situa-se 50% acima do valor de referência; 

iv. os parâmetros óptimos parecem assumir valores que dão ao conjunto uma con- 
figuração com características muito especiais, como a seguir se perceberá. 

Repare-se que os parâmetros óptimos parecem aproximar-se dos valores h = \f2 = 1.41421 . . . 
e /3 = 0. Coloca-se então a seguinte questão: 

JVão serão esses os valores exactos dos parâmetros óptimos? 

A representação gráfica da função R{h, (3) através de curvas de nível. Figura [2] está em 
perfeita concordância com essa possibilidade — repare-se que as curvas de nível parecem 
perfeitamente centradas no ponto de coordenadas (\/2, 0), assinalado na figura por um 
Veja-se também, na Figura [sj o gráfico da resistência R{h) para /? = O, onde é igualmente 




1.42 



(a) (b) 
Figura 2: Curvas de nível da função R{h,(3). 



perceptível uma surpreendente elevação da resistência quando h ^/2. Calculou-se então 
numericamente a resistência da cavidade Q^'»'' com os valores exactos h = y/2 e /? = O, tendo 
o resultado confirmado o valor 1.49650. 

Existe ainda uma razão adicional que sugere também uma resposta afirmativa à questão 
formulada. A forma do conjunto íí^^íi./s ^gm h = V2 e p = O é um caso particular a que 
estão associadas características especiais que poderão justificar o elevado valor de resistência 
que apresenta. As duas secções da forma são arcos equivalentes de duas parábolas de eixos 
horizontais e concavidades voltadas uma para a outra — ver Figura |4j Mas a particularidade 
da configuração reside no facto do eixo das parábolas coincidir com a linha de entrada da 
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Figura 3: Gráfico da resistência R{h) para /? = 0. 



cavidade (eixo dos x), e o foco de cada uma situar-se no vértice da outra. Repare-se que, sendo 
a equação da parábola do lado direito dada por x = Q^qÍu), obtém-se = — 4(x — 1/2), 
confirmando que o vértice é em (1/2,0) e que o foco dista deste uma unidade. Desse modo, 
por razões de simetria, os focos e os vértices das duas parábolas situam-se nos extremos da 
abertura da cavidade, (—1/2,0) e (1/2,0), tal como se ilustra no esquema da Figura|4]3. 




(a) (b) 
Figura 4: Forma 2D (quase) óptima — futuramente conhecida por Dupla Parábola. 



O elevado valor da resistência desta forma de cavidade está ainda expresso, embora de 
modo implícito, no gráfico da função integranda da funcional apresentado na Figura |5j 
Trata-se de uma superfície com uma forma bastante suave em quase todo o seu domínio e, 
mais importante, quase não varia com x e parece descrever a forma da função 2 cos ip na 
direcção do eixo dos ip. Estes atributos, quando verificados na sua plenitude, estão associados 
ao maior valor de resistência admissível, R = 1.5j^ 

^Se a função integranda tomar a forma G(x,(p) = 2cos<^, de (|3| obtém-se R — 1.5. 
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Figura 5: Gráfico da função integranda G{x, íp) = {1 + cos (99"^ (x, íp) — íp)) cos ip. 



Esta forma de cavidade parece tratar-se efectivamente de um caso muito particular. Con- 
trariamente ao que se passou com todas as outras formas anteriormente estudadas, a função 
integranda aparenta uma forma bastante suave, apresentando apenas pequenas irregulari- 
dades para ângulos ip de pequena amplitude. Atendendo a essa característica, e tendo em 
conta que G{x, ip) quase não depende de x, a resistência foi calculada, para esta forma em 
particular, usando a regra de Simpson 1/3 na integração em ordem a ip. A dupla integração 
na equação ([s]) foi então aproximada numericamente pela seguinte expressão: 

R=-AxAíp ^ Wk {1 + cos {íp^{xi,(pk) - ípk)) cos ípk, (9) 

j=Ar^/2+l k=l 

com Wk = 2 para k ímpar e Wk = 1 para k par, Xi = —l/2 + {i — 1/2) Ax, Ax = 1/Nx, 
ipk = — vr/2 -|- kA^p e Api = Ti/Nip. e A'^^ são o número de subintervalos a considerar 
na integração das variáveis x e, ip (ambos números pares), respectivamente, e Ax e Aip os 
incrementos para as correspondentes variáveis discretas. Dado que a forma íí^v^-o apresenta 
simetria horizontal, o primeiro somatório da expressão considera apenas a segunda metade 
do intervalo de integração da variável x. 

Para que passe a ser facilmente referenciável, esta forma de cavidade (Figura [4^) passará, 
a partir de agora, a ser designada simplesmente por ^^Dupla Parábola" . Assim, no contexto 
deste trabalho, o termo "Dupla Parábola" deve ser sempre entendido como o nome da cavidade 
cuja forma é descrita por duas parábolas que, para além de serem geometricamente iguais, 
encontram-se "encaixadas" na disposição peculiar que referimos. 

Uma vez que a resistência da Dupla Parábola assume um valor já muito próximo do seu 
limite teórico, numa derradeira tentativa de se vir a conseguir atingir esse limite, resolveu-se 
estender ainda o estudo a outras classes de funções g{y) que admitissem a Dupla Parábola 
como caso particular ou que permitissem configurações próximas dessa forma quase óptima. 
Em todos esses casos os melhores resultados foram invariavelmente obtidos quando a forma 
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das curvas se aproximou da forma da Dupla Parábola, sem nunca terem superado o valor 
R = 1.4965. Começou-se por considerar funções g{y) seccionalmente quadráticas, incluindo 
curvas splines, sem que se conseguissem resultados interessantes; apenas para funções g{y) de 
2 ou 3 segmentos foi possível aproximarmo-nos da resistência e da forma da Dupla Parábola. 
Consideraram-se também funções g{y) cúbicas e biquadrática^ mas em ambos os casos o 
processo de optimização aproximou-as a curvas de ordem quadrática, com os coeficientes de 
maior ordem a tomarem valores quase nulos. Estudou-se o problema na classe das secções 
cónicas, considerando-se, para faces laterais da cavidade, dois arcos simétricos quer de uma 
elipse quer de uma hipérbole. Também nestes casos os arcos assumiram uma forma muito 
próxima dos arcos de parábolas. 

Sendo a Dupla Parábola a melhor forma encontrada, e tratando-se de uma forma quase 
óptima, na secção que se segue é objecto de um estudo aprofundado, de natureza essencial- 
mente analítica, onde se tentam perceber as razões do seu bom desempenho. 

4 Caracterização das reflexões na forma "Dupla Parábola" 

Cada uma das ilustrações da Figura |6] reproduz, para a "Dupla Parábola", uma trajectória 
concreta, obtida com o nosso modelo computacional. E reconfortante verificar que, à excepção 




(a) X = 0.45, if = 75°. (b) x = 0.45, 99 = 55°. (c) x = 0.45, ip = 35°. 




(d) x = 0.3, if = 75°. (e) x = 0.0, ip = 35°. (f) x = 0.48, tp = 5°. 
Figura 6: Exemplo de trajectórias obtidas com o modelo computacional. 

^Nas curvas biquadráticas, o ponto de intercepção da trajectória da partícula com a fronteira da cavidade é 
calculado resolvendo uma equação do 4° grau. As raízes dessa equação foram obtidas numericamente usando 
o método descrito em [8]. As equações de ordem inferior foram sempre resolvidas recorrendo às conhecidas 
fórmulas resolventes. 
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da última trajectória, em todas as restantes a partícula sm'ge à saída da cavidade com uma 
velocidade quase invertida relativamente àquela que foi a sua velocidade de entrada. Este é 
o "sintoma" que caracteriza inequivocamente uma cavidade de óptimo desempenho. Mesmo 
no caso da trajectória da ilustração (f), a direcção da velocidade de saída parece não andar 
muito distante da de entrada. 

Se analisarmos as cinco primeiras ilustrações, verificamos existir algo em comum no com- 
portamento da partícula: para descrever a trajectória, a partícula é sempre sujeita a três 
reflexões. Esta parece ser uma característica determinante para a aproximação dos ângulos 
de entrada e de saída. Se, por exemplo, imaginarmos três trajectórias com configurações 
próximas, respectivamente, das trajectórias (a), (b) e (c), mas com a diferença de não 
possuírem a terceira reflexão, o resultado será completamente diferente, como facilmente se de- 
preende das ilustrações. Embora esta convicção seja por enquanto de natureza essencialmente 
empírica, os resultados do estudo que se segue vão no sentido de confirmar que uma parte 
muito significativa das trajectórias "benignas" — aquelas em que os vectores velocidade de 
entrada e de saída são quase paralelos; chamemos-lhes assim por representarem contribuições 
positivas na maximização da resistência — subentendem exactamente três reflexões. 

Vamos agora tentar interpretar outro tipo de resultados obtidos com o nosso modelo com- 
putacional, começando pela representação gráfica do valor da diferença entre os ângulos de 
entrada e de saída — (p~^{x, ip)), em função da posição de entrada a; e do ângulo de entrada 
(/9 — ver Figura [7j O gráfico foi produzido usando 100 valores, uniformemente distribuídos, 
para cada uma das variáveis x S (0,0.5) e G (—90°, 90°) — dada a simetria da cavidade 
em relação ao eixo x = O, obter-se-ia um resultado equivalente para o intervalo x G (—0.5,0). 
Este gráfico mostra-nos que, independentemente do valor da posição de entrada x, a diferença 




Figura 7: Representação gráfica da diferença angular ip — 99"'" (x, 99). 

angular é aproximadamente nula para ângulos de entrada 99 de amplitude elevada, apresen- 
tando valores menos interessantes quando íp se aproxima de zero. Portanto, começa-se a 
perceber que as trajectórias "benignas" têm origem essencialmente em ângulos de entrada de 
amplitude elevada. 

Os gráficos das figuras |8]e|9] que se seguem no nosso estudo, foram produzidos com 10.000 
pares de valores (x, 99), gerados por um processo aleatório de distribuição uniforme. O primeiro 
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se nas proximidades da diagonal (p = 99^, o que é revelador de um bom comportamento 
da cavidade. Também com estes resultados se confirma que a resposta da cavidade se vai 
deteriorando quando ip se aproxima de zero. 

Vejamos o segundo gráfico, expresso na Figura [9] onde se encontra representada no plano 
cartesiano a distribuição dos pares {(p,ymax), sendo ymax a altura atingida pela partícula. 
Pela análise da figura, constatam-se essencialmente duas coisas: que a altura atingida pela 
partícula quase não depende da posição de entrada x, dependendo essencialmente do valor 
do ângulo de entrada íp, em especial para ângulos elevados — esta relação é perceptível na 
Figura |6] quando comparamos a ilustração (a) com a (d) e a (c) com a (e) — ; e que a altura 
que a partícula pode atingir aumenta com a diminuição da amplitude do ângulo de entrada. 

Comparando agora os gráficos [8] e |9] percebe-se que existe uma forte correlação entre o 
desempenho da cavidade e a altura atingida pela partícula: quanto mais elevada é a altura 
atingida, maior é a dispersão dos valores em relação à diagonal íp = p^. Este facto leva- nos 
a supor que a forma do contorno da cavidade é a ideal junta da sua base, mas vai deixando 
de o ser à medida que nos aproximamos do topo da cavidade. 

A zona da cavidade que interage com a partícula parece não ser o único factor a condi- 
cionar a similitude entre os ângulos íp e ip~^. Se repararmos no gráfico [sj parece existir uma 
perturbação adicional no comportamento da cavidade quando a amplitude do ângulo de en- 
trada é inferior a cerca de 20°, que leva a que alguns pares {p,p^) fiquem, em relação aos 
restantes, mais dispersos e mais distantes da diagonal p^ = ip. Chamámos já a atenção para 
a possível importância das 3 reflexões no grau de aproximação verificado nos ângulos (p e p^. 
Ocorre-nos por isso a seguinte questão: não será precisamente o número de reflexões que, ao 
se diferenciar das 3 ocorrências, interfere tão negativamente no comportamento da cavidade? 
As investigações que se seguem vão demonstrar, entre outras coisas, que esta nossa suspeita 
tem fundamento. 
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A confirmar-se a nossa última conjectura, e atendendo ao aspecto da distribuição dos pares 
no gráfico [sj é expectável que exista para (/? um intervalo {—(po,(po) C (—20°, 20°) 
que abarque todas as trajectórias que não sejam formadas por 3 reflexões. A existência e 



identificação desse intervalo são garantidas com a demonstração apresentada na subsecção 4.1 
do postulado que se segue: 

Para ângulos de entrada Lp superiores (em valor absoluto) a (po = 
arctan (^^^ ~ 19.47°, o número de reflexões a que a partícula é sujeita no interior 
da cavidade é sempre igual a três, e ocorrem alternadamente nas faces esquerda e 
direita da cavidade, para qualquer que seja a posição de entrada. 

De forma a verificarmos que as deduções que fizemos estão efectivamente em concordância 
com os resultados numéricos do modelo computacional desenvolvido, apresentamos mais dois 



gráficos. Figuras 10 e 11, ambos produzidos com 10.000 pares de valores {x,(p), gerados 
aleatoriamente com distribuição uniforme. Como pode ser observado na Figura [TÕ] todas 
as trajectórias com 4 ou mais reflexões (colisões), entre as 10.000 consideradas, aconteceram 
dentro do intervalo {—ípo, (fo)- Fora desse intervalo (para \(p\ > (po) as trajectórias são sempre 
de três reflexões. Adicionalmente, podemos também verificar não existir qualquer trajectória 
com menos de três reflexões. Provaremos na subsecção |4 . 2 1 ser esta também uma característica 
da cavidade. 



Não menos elucidativo é o resultado apresentado no gráfico da Figura 11 onde difer- 
enciámos dos restantes os pares (ip, ip^) que estão associados a trajectórias com 3 reflexões. 
Constata-se, tal como suspeitávamos, que os pares (9?, ip~^) que se encontram mais dispersos e 
mais distantes da diagonal ip^ = ip, estão todos eles associados a trajectórias com 4 ou mais 
reflexões. 

Das conclusões a que chegámos podemos de imediato retirar o seguinte corolário: mesmo 
nas condições que conduzem ao pior desempenho da cavidade — 4 ou mais reflexões — , a 




Figura 9: Distribuição no plano cartesiano dos pares {íp,ymax), sendo Umax a altura atingida 
pela partícula. 
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diferença angular \ip — 93"*" |, por maior que seja, nunca será superior a 2ipQ ~ 38.94°, valor 
que é bastante inferior ao maior ângulo que é possível formar entre dois vectores (180°). 
A demonstração deste corolário é simples: como uma trajectória de 4 ou mais reflexões 
está sempre associada a um ângulo de entrada —ipQ < ip < ipo, o ângulo de saída situar- 
se-á necessariamente no mesmo intervalo; tendo em conta a propriedade de reversibilidade 

I 



I 
I 



1^ "fio 



-80 -60 -40 -20 O 20 40 60 80 
Figura 10: Distribuição no plano cartesiano dos pares {(p,nc), sendo nc o n° de colisões. 




-20 -15 -10 -5 o 5 10 15 20 



Figura 11: Distribuição dos pares (íp,(p~^): (•) pares associados a trajectórias com 3 reflexões; 
(o) pares associados a trajectórias com 4 ou mais reflexões. 
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associada à lei de reflexão que rege as reflexões, se por absurdo admitíssemos \ip^\ > (po, ao 
invertermos o sentido de marcha da partícula, passaríamos a ter uma trajectória de mais de 
3 reflexões com um ângulo de entrada ip~^ situado fora do intervalo {—(fio, ^o), o que entraria 
em contradição com o postulado inicial. Os resultados apresentados na Figura 11 estão em 
concordância com este novo facto, uma vez que todos os pares (^^) com mais de três 
reflexões parecem situar-se no interior da região (—990,920) x {—^o^V^o)- 
Resumindo: 



• A altura que a partícula atinge no interior da cavidade interfere negativamente na 
aproximação dos ângulos ip e </?"'", levando-nos a pensar que a forma da cavidade é a 
ideal, ou muito próxima da ideal, junto à sua base, mas vai perdendo qualidades em 
zonas mais elevadas da cavidade; 

• O maior desfazamento entre íp e ocorre em trajectórias de 4 ou mais reflexões; 

• Verifica-se, em todo o intervalo de variação de íp, uma grande predominância das tra- 
jectórias com 3 reflexões; 

• Não existem trajectórias com menos de 3 reflexões; 

• O ângulo crítico ipo tem o valor: ipç, = arctan (^^^ — 19.47°; 

• Fora do intervalo {—ipq,ípq), todas as trajectórias são de 3 reflexões; 

• Em trajectórias com 4 ou mais reflexões, a diferença angular é delimitada por 2(po: 

\íp - íp+\ < 2lpq. 



4.1 Uma condição suficiente para a ocorrência de três refiexões 
Teorema 1. Para ângulos de entrada (p superiores (em valor absoluto) a ipo = arctan 



19.47°, o número de reflexões a que a partícula é sujeita no interior da cavidade Dupla 
Parábola é sempre igual a três, e ocorrem alternadamente nas faces esquerda e direita da 
cavidade, para qualquer que seja a posição de entrada. 

Para demonstrarmos o teorema que acabámos de enunciar, estudaremos a trajectória de 
uma partícula passo a passo, desde o momento em que entra na cavidade até ao momento 



em que a abandona, e usaremos as ilustrações das figuras [T2j|T4] para nos auxiliarem nessa 
demonstração. 

Considere-se uma partícula que entra na cavidade em (x,0), com o vector velocidade a 
formar um ângulo (p com o eixo vertical, tal como se encontra representado nas ilustrações da 



Figura 12 , onde assumimos que o eixo de simetria da cavidade é o eixo dos y e que a sua base 



AqAi assenta no eixo dos x. Assim, a posição da partícula à entrada da cavidade assume 
apenas valores no intervalo (—5, 5) x {0}. 

Dada a simetria da cavidade em relação ao seu eixo vertical, será suficiente analisar o seu 
comportamento para (po < ip < 90°. As conclusões a que chegarmos serão assim igualmente 
válidas para —90° < ip < —ípo- 

Analisemos então em detalhe e separadamente cada um dos subtrajectos que compõem 
toda a trajectória descrita pelo movimento da partícula no interior da cavidade. 
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Subtrajecto BqB\ 

Para (p > (Çq, temos a garantia de que a primeira reflexão ocorre na curva parabólica 
do lado esquerdo da cavidade, tal como pode ser facilmente deduzido a partir da ilustração 
(a). Para que a partícula colida na curva esquerda bastaria que o ângulo fosse superior 

a arctan^^^^, grandeza que tem como majorante içq = arctan^^^. Teremos assim o 
trajecto inicial da cavidade representado na ilustração (a) pelo vector BqBi. 



Subtrajecto B1B2 
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Após colidir em Bi, de acordo com a lei de reflexão, a partícula segue pelo trajecto B1B2. 
> 

Demonstremos que B1B2 tem o sentido ascendente — ilustração (a). Tracemos o segmento de 

recta A1A2, paralelo ao trajecto inicial da partícula BqBi e que passe pelo foco da parábola 

esquerda {Ai). Pela propriedade focal dessa parábola, uma partícula que tome o subtrajecto 

A1A2, após a reflexão em A2, seguirá numa direcção horizontal A2A3 (prosseguindo depois o 

seu trajecto, após nova reflexão, em direcção ao foco Aq da segunda parábola). Ocorrendo a 

primeira reflexão da partícula em Bi , um ponto da curva necessariamente posicionado abaixo 
> 

de A2, o trajecto B1B2, que seguirá de imediato, será no sentido ascendente, pois a derivada 
^ da curva nesse ponto (Bi) é superior à derivada em A2, onde a trajectória que se seguia 
era horizontal. 

> 

Embora saibamos já que -B1-B2 tem o sentido ascendente, ainda nada nos garante que 
a segunda reflexão aconteça necessariamente na parábola do lado direito. Se conseguirmos 
veriflcar que para = ipo a segunda reflexão é sempre no lado direito, para qualquer que 
seja a posição de entrada x, então, por maioria de razão, o mesmo sucederá para qualquer 
valor (p > ipQ. Esta premissa pode ser facilmente aceite com o auxflio da ilustração (a) da 
Figura 13 para qualquer valor de (p > ipo, com a primeira reflexão num dado ponto Bi, é 
sempre possível traçarmos uma trajectória para ip = ipo que apresente a primeira reflexão 
no mesmo ponto Bi; sendo a segunda reflexão na curva do lado direito para o caso íp = ipo, 
necessariamente o mesmo acontecerá para a trajectória com (p > ipo, pois o ângulo de reflexão 
será menor neste segundo caso, tal como se ilustra na flgura. Por conseguinte, bastar-nos-há 
provar para ip = (po, que a segunda reflexão ocorre sempre na parábola do lado direito, para 
que o mesmo flque provado para qualquer que seja o ip > ipo. 





(a) (b) 

Figura 13: Ilustrações para estudo da segunda reflexão. 



Na ilustração (b) da Figura 13 encontra-se representada a trajectória até à segunda re- 
flexão de uma partícula com ângulo de entrada ipo {BqBi e B1B2). Como se depreende da 
ilustração, a reflexão B2 só acontecerá na curva do lado esquerdo se o ângulo a for menor que 
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a' . Determinemos o valor dos dois ângulos. 

Sendo (xi,yi) as coordenadas do ponto teremos tan(Qí') = ; logo 




a' = arctan j^ -(W4-l/2) ^ ^ ^^^^^^ f \2-y^)/A ~\^ ^ ^^^^^^ j j _ ^^^^ 

Para chegarmos ao valor de a resolvemos o sistema de três equações, de incógnitas a, O e P, 
que se retira directamente da geometria da própria figura 

a + p + 9 = TT 

arctan (|yi) + (po + 9 = ^ 

A linha tangente à curva em Bi faz com a vertical um ângulo cuja tangente tem por valor 
a derivada g| da curva nesse ponto (em y = yi), onde g| = ^(3:2/^ — 5) = ^y- Por isso, esse 
ângulo surge representado na terceira das equações pela grandeza arctan (^yi)- 

Resolvendo o sistema, obtém-se para a o seguinte resultado 

a = ipQ + 2 arctan = arctan ( ) ~'~ ^ arctan • (H) 



Provemos finalmente que a > a' , para qualquer que seja yi G (O, v2)- Das equações (11 ) 



e (10), será equivalente a provarmos 



arctan [ 1+2 arctan f — 1 > arctan í ^ ] . 
[4) \2) \ ^ ) 

Dado que O < yi < \/2, ambos os membros da inequação representam ângulos situados no 
primeiro quadrante do círculo trigonométrico. Por isso podemos manter a inequação para 
a tangente dos respectivos ângulos. Aplicando a tangente a ambos os membros, depois de 
efectuadas algumas simplificações trigonométricas, chegamos á seguinte relação 

1 4^/2 - ^/2yl + 16yi ^2 + yi 
4 4-y2_y2yi ^ 4 ' 

que, com simplificações algébricas adicionais, toma a forma 

yi {V2yi + 14 + yf) 



y? - V2yi 



> 0. 



Como O < yi < \/2, facilmente se constata que tanto o numerador como o denominador da 
fracção presente nesta última inequação são grandezas positivas. Logo a > a' , que contraria a 
condição que era necessária para que a reflexão B2 ocorresse na curva do lado esquerdo, ficando 
assim provado, como pretendíamos, que em nenhuma situação a reflexão B2 da ilustração (b) 
da Figura [13] acontece na curva do lado esquerdo. Por maioria de razão, podemos então 
também concluir que o mesmo sucede para qualquer que seja tp > ipo- a segunda reflexão da 
partícula ocorre sempre na parábola do lado direito. 
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Subtrajecto B2B3 



Demonstremos que o subtrajecto -B2-B3 tem o sentido descendente — ilustração (b) da 



Figura 12 Imagine-se, para o efeito, um subtrajecto AqA2, paralelo a B1B2 e que passe no 



foco Ao. O subtrajecto A2A3 que se seguiria à reflexão em A2 — um ponto da parábola 
do lado direito situado abaixo de B2 — seria horizontal. Sendo a derivada da curva em A2 
superior ao valor da derivada em B2, o subtrajecto B2B3 será necessariamente de natureza 
descendente. 

Ainda que já saibamos que o subtrajecto é descendente, ainda não mostrámos que esse 
subtrajecto em nenhuma situação conduz a partícula directamente para a saída da cavidade. 
Segue-se então a demonstração de que a reflexão -B3 ocorre sempre numa posição superior 



a Ao — ilustração (c) da Figura 12, Tracemos A2B2, um segmento de recta horizontal que 
passe no ponto de reflexão B2. Se a partícula seguisse esse trajecto, colidiria no mesmo ponto 
B2, mas dirigir-se-ia para Aq. Logo, pela lei de reflexão, B3 terá que estar acima de ^o, pois 
B1B2 faz um ângulo com o vector normal à curva em B2 menor que o formado pelo segmento 



Subtrajecto B3B4 

Vamos agora mostrar que o subtrajecto que se segue à reflexão em B^ cruza o segmento 

Tracemos 
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^0^1 1 isto é, direcciona-se para fora da cavidade — ilustração (d) da Figura 

então A2B3, um segmento de recta horizontal que passe no ponto de reflexão B3. Se a partícula 

seguisse esse trajecto, colidia em e dirigir-se-ia para Ai. Logo, pela lei de reflexão, a recta 

> ^ 

onde assenta o subtrajecto B3B4 terá necessariamente que passar abaixo de Ai, pois B2B3 

faz um ângulo com o vector normal à curva em B^ maior que o formado pelo segmento A2B3. 

Mostrámos que o subtrajecto cruza o eixo dos x num ponto situado à esquerda de Ai, mas 

ainda não mostrámos que ocorre à direita de Aq. Para tal, teremos que demonstrar que a 

terceira é a última das reflexões, isto é, que em nenhuma situação ocorre uma quarta reflexão 

na parábola do lado esquerdo. Segue-se essa demonstração, de todas a mais complexa e a 

mais demorada. 

Para provarmos que a seguir à terceira reflexão não ocorre qualquer outra colisão na 
parábola esquerda, vamos mostrar que uma quarta colisão — representada por B4 na ilus- 



tração (a) da Figura 14 — tem sempre origem num ângulo de entrada (f inferior a (fQ. Iremos 



assim estudar o trajecto da partícula na ordem inversa à sua progressão: começamos por 
admitir a existência do subtrajecto -B3-B4 da ilustração (a) e analisaremos as suas implicações 
em todo o trajecto precedente. 



Na ilustração (a) da Figura 14 encontram-se representados os subtrajectos B2B3 e -63-64. 



Comecemos por relacionar 02 com «3, os ângulos que os vectores -B2-B3 e B^Bí^, respectiva- 
mente, formam com o eixo vertical. Para o efeito resolvemos o sistema de três equações, de 
incógnitas «21 ^3 e /32, que se retira da geometria da figuraj^ 

6-i = 1^2 + «3 
arctan (^^3) /Ja = f 
«2 + ^3 + /?2 = vr 

obtendo-se 

«2 = 2 arctan — «3; 



''As variáveis denotadas por Xi e yi, com i — 1,...,4, representam as coordenadas do i-ésimo ponto de 
reflexão, identificado por Bi. 



19 



em que arctan (52/3) é o ângulo que a recta tangente à curva em faz com a vertical 
o declive da recta tangente é dado por ^ \y^y^ = - 5) \y=y^ 

ângulo as pode ser expresso da seguinte forma 



«3 = arctan 



X3 - X4 



arctan 



1„2 _ 1,,2 
4^3 4^4 



arctan 



= 2^3- Por sua vez, o 
y3 + 2/4 



.v^.-Vij \ y3-y4 j V 4 

o que nos permite escrever «2 em função unicamente das ordenadas 2/3 e 2/4 dos extremos do 
> 

vector -B3-B4, 



o 4- /'ysA , / 2/3 + 2/4 

«2 = 2 arctan ( — | — arctan 



, , . (12) 

Para conseguirmos provar o que pretendemos — impossibilidade de ocorrência da reflexão 
S4 — precisamos de encontrar um minorante para a ordenada da posição onde ocorre cada 
uma das quatro reflexões, ou seja, determinar {2/*, 2/2; 2/3; 2/4}) tal que 



2/1 > 2/í, 2/2 > 2/2> ys > yh Vi ^ Vl^ V(y?,x) G (^930, 1) X ^- 



1 1 
2' 2 



(13) 



Facilmente se percebe que 2/4 = 0. Vamos então determinar os outros três minorantes, 
começando por 2/2- 



Sabemos que O < 2/4 < 2/3! logo, de (12) retiramos que 
arctan < 02 < 2 arctan 



arctan 



2/3 



Atendendo a que «2 se situa no primeiro quadrante do círculo trigonométrico, podemos man- 
ter as desigualdades para a tangente dos respectivos ângulos. Após algumas simplificações 
algébricas, obtemos 

ys ^ ^ 2/3(12 + 2/1) 



2 < tan(Q2) < 



(14) 





(a) (b) 

Figura 14: Ilustrações para estudo de uma hipotética quarta reflexão. 
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A equação da recta que liga B2 a B3 toma a forma 

com m = tan(a2) 6X3 = — ^. Como estamos interessados em encontrar a ordenada do 
ponto de intercepção dessa recta com a curva parabólica situada no lado direito, de equação 

1 2 1 
^ = -4^ +2' 

temos de resolver a equação de segundo grau, na variável y, que resulta da eliminação da 
variável x por combinação das duas equações anteriores. A ordenada da segunda reflexão, 
sendo a raiz positiva da equação, toma a forma 



y2 = — 2m + Y 4m^ + ^my^ — y| + 4 . 

A grandeza y2 é expressa em função de duas variáveis, m e y^, que como sabemos assumem 
apenas valores positivos. De forma a aceitarmos mais facilmente as deduções que iremos fazer 
no encalço de y2, imaginemos, sem qualquer perda de generalidade, que 2/3 é um valor fixo. 
Comecemos por mostrar que a derivada de 2/2 em ordem à variável m. 



(15) 



dy2 2 (^2m + 2/3 - \/4m2 + 4my3 - y| + 4^ 
dm y/4;^^2Tp4;^^^J^rr^^|Tp4 

tem um valor negativo para qualquer que seja o valor de ^3. Como 7/3 < V2, forçosamente 
2/3 < 4, logo os dois radicandos {Am? + 4my3 — 2/3+4) presentes na equação ( 15 ) têm sempre 
um valor positivo. A restrição 2/3 < \/2 permite-nos ainda deduzir sucessivamente as seguintes 
desigualdades 

2/i < 2 4^ 22/i < 4 4^ 2/1 < 4 - 2/i <^ ^m^ + ^my^ + 2/| < 4771^ + Amy^ + 4-2/1 



<^ {2m + 2/3) < 4m + 4mys - 2/3 + 4 44> 2m + 2/3 < y 4m2 + 4m2/3 - j/f + 4 . 

Esta última desigualdade confirma que ^ < O, para qualquer que seja o 2/3. Assim, o valor 2/2 

é tanto menor quanto maior for o valor de m. Como se depreende de ( 14 ), m < M = ^^Í^?^^^^ 
logo 



2/2 > -2M + ^471^2 + 4M2/3 - 2/1 + 4 . 
Substituindo M, obtém-se, após algumas simplificações, 

3 1,1 



2/2 > /(ys), com /(2/3) = -^ys - + g V^^^^i + 402/1 + 2/^ + 256 . (16) 

Para encontrarmos o valor mínimo de 7(2/3) começamos por derivar, 

2722/3 + 8O2/Í + 32/i - (12 + Syj) ^212yl + 402/| + 2/f + 256 
d2/3 8^/2722/1 +4O2/I + 2/Í + 256 

Sendo os radicandos claramente positivos, apenas temos que nos preocupar com o numerador 
da fracção. Encontrarmos as raízes da função derivada 3^/(2/3) equivale por isso a resolvermos 
a equação 

(2722/3 + 8O2/Í + 32/i)' = (12 + 32/Í)' (2722/i + 402/| + 2/3 + 256) , 
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que pode ser simplificada na seguinte: 

2304 - 1024y| - 992y| - 160y^ - 3yl = 0. 
Esta equação polinomial tem uma única raiz real positiva, de valor 

m = ^^-51 + 6^79 ~ 1.017, 

significando que /(ys) tem um mínimo global em y^, pois, como mostramos a seguir, 
^fius) > O e a função não apresenta outros pontos de estacionaridade. 

Mostremos então que ^/(ya) > O, com 

d^ ^ _ 18240y| + 2960yf + 180yf + 3y;^° + 34816 + 30720y| 
dyl 4(272y| + 40y| + yf + 256) § 

(816yi + 120y| + 3y| + 768y3) V272y| + 40y| + yf + 256 



4(272y| + 40y| + yg + 256) § 



(17) 



Mostrarmos que x^/í^s) > O equivale a mostrarmos que o numerador da primeira fracção é 



superior ao numerador da segunda, em (17). Depois de elevarmos os dois termos ao quadrado, 
chegamos à inequação 

-3y;^^ + 36y;^'^ + 9464y;^2 + 191616y;^° + 1514112y^ + 5817344yf 
+13535232y^ + 15532032y| + 9469952 > 0. 

Facilmente comprovamos a veracidade desta relação, dado que temos um único termo negativo 
(— 3y3^) que, por exemplo, é inferior em valor absoluto ao termo constante (9469952), 

y^<^/2^ 3yl^ < 768 < 9469952. 

Fica assim completa a demonstração de que /(ys) tem um mínimo global em ys, de valor 



fim) = -V-51 + 6V79. 



Portanto, de (16), concluímos finalmente que 



2/2 >y2 = gV + 6V79~ 1.356. (18) 
Está então encontrado um minorante para a altura da segunda reflexão B2 (ilustração (a) da 



Figura 14). Determinemos agora y^, um minorante para a altura da terceira reflexão -B3. 

Para que a reflexão B2 ocorra na parábola do lado direito é necessário que o ângulo «2 
seja maior que o ângulo formado entre o eixo vertical e o segmento de recta que une com 
o vértice superior da cavidade. 
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Esta inequação, conjuntamente com a segunda relação de desigualdade de ( |14[ ), permite- nos 

< ^ ^ < ' 

de que resulta a inequação 

Vi 



+ 8y3 + > 0. 



Como o polinómio da inequação tem derivada positiva e admite uma única raiz real, con- 
cluímos de imediato que a mesma constituiu um limite inferior para ys, sendo esse limite 

1 / \ ''^ 8 
y* = - (54\/2 + 6\/54f))' ^- 0.670. (20) 

(54^/2 + 6^/546) ^ 



3 

Resta-nos determinar yj, um minorante para o valor de yi — ordenada onde ocorre a 



primeira reflexão. Para o efeito recorremos à ilustração (b) da Figura 14 que nos dá uma 
representação mais pormenorizada da parte da cavidade onde ocorrem as duas primeiras 
reflexões, Bi e B2. O esquema apresentado foi construído contando que a primeira reflexão 
(Bi) ocorre num ponto mais elevado do que o da terceira reflexão (-B3). E de facto essa a 
situação. Isso mesmo pode ser comprovado mostrando que 02 é sempre menor que o ângulo 
formado entre o vector normal à curva em B2 e o eixo vertical, ou seja 

vr /l 
«2 < 2 ~ arctan ( -y2 



Pegando, em (14), no limite superior de tan(a2) e tendo presente que y < \/2, construímos a 



seguinte sequência de desigualdades que comprova o que se pretende: 

~51.06° ~54.74° 

/y3(l2 + yi)\ ^ (7V2\ tt Í V2\ tt íl 

a2 < arctan I ^-^ 1 < arctan I ] 2 ~ l ] ~2~ l 

Tentemos agora encontrar y\. Podemos deflnir y\ como sendo a ordenada do ponto de 
intercepção da parábola esquerda com a semi-recta de origem no ponto B2, posicionado o 
mais abaixo possível (2/2 = 2/2)1 ^ declive igual ao maior valor permitido para o declive 
da trajectória que antecede B2 {B1B2). A equação da recta que liga Bi a B2 toma a forma 

X = m{y - 2/2) + X2, 

com m = tan(ai) e X2 = — jí/I + 2- Como estamos interessados em encontrar o ponto de 
intersecção dessa recta com a curva parabólica situada no lado esquerdo, de equação 

1 2 1 

X = -y , 

4^ 2 ' 

teremos que resolver a equação de segundo grau, na variável y, que resulta da eliminação 
da variável x por combinação das duas equações anteriores. Embora tenha duas raízes reais 
positivas, apenas estamos interessados na menor delas, que toma a forma 



yi = 2m — J Am? — 4my2 — 2/1 + ^ • (^1) 
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Como dissemos, se fizermos y2 = Hq, ^ impusermos o declive máximo para a recta, que 
nas equações anteriores equivale a considerarmos m mínimo, obtemos yi = y*- Sendo m = 
tan(Q;i), devemos determinar o valor de ai através do sistema de equações 

«1 = 202 + 02 

arctan (51/2) + «2 + é'2 = f 



que se retira da geometria da ilustração (b) da Figura 14 Obtem-se 



oi = vr — 2 arctan 



«2- 



(22) 



Desta última igualdade, de (14), e dado que y2<v2ey3<v2, deduzimos que 

'y3(l2 + yi)' 



ai > TT — 2 arctan 



y2 



arctan 



16 



> TT — 2 arctan 



V2' 



arctan 



7 



logo. 



m = tan(ai) > tan [ vr — 2 arctan 



V2' 



arctan 



7 
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V2. 



Se em (21 ) fizermos m = |i\/2 e 7/2 = 2/2 obtemos então o limite inferior para yi 



10 



90 



yl = ^V2- V 444498 -33120\/2V -51 + 6V79 -38400\/79 ~ 1.274. 



(23) 



Resumindo, 



iyl,y*2,ylyl) ^ (1-274, 1.356, 0.670, 0). 



Com o auxílio da ilustração (b) da Figura 14 vamos, por fim, analisar o ângulo de entrada 
íf da partícula. Com o sistema de equações (^yi = g| \y=yi) 

29i + í/9 + ai = vr 
arctan {\yi) + + 6*1 = | 



e com as igualdades (22) e (12) obtemos, sucessivamente, 

'yv 



í/3 = «1 — 2 arctan 
í/3 = vr — 2 arctan 



2 ; ' 



yi 



yi 



2 arctan 



í/9 = vr — 2 arctan I — ) ~ 2 arctan l — ) ~ 2 arctan l — ) + arctan 



2/2 



y2 



«2, 



(24) 



Pegando em (24), de (19) deduzimos que 

'2/1 



í/? < vr — 2 arctan 



2 arctan 



2/2 



arctan 



V2 + 



2/3 



De acordo com as definições (13) e com os valores encontrados em (18), (20) e (23), podemos 
concluir que 



í/3 < vr — 2 arctan 



2/1 



2 arctan 



arctan 



V2 + 2/3 



19.18°, 
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ou seja 



(/:)< V?o - 19.47°. 



Com isto podemos finalmente concluir que é impossível termos uma quarta reflexão, pois, 
a acontecer, a partícula teria que ter entrado na cavidade com um ângulo 99 < (/jq, como 
acabámos de mostrar — algo que contraria a nossa imposição inicial, ip > ipQ. Como a 
cavidade apresenta simetria em relação ao seu eixo vertical central, a conclusão a que chegamos 
é igualmente válida para íp < —(fo, ficando assim provado o que pretendíamos (Teorema [T|: 

Para \(p\ > (po, ocorrem sempre três refíexões, alternadamente nas faces es- 
querda e direita da cavidade Dupla Parábola. 

4.2 Número mínimo de reflexões 

Teorema 2. Qualquer partícula que entre na cavidade Dupla Parábola descreve uma tra- 
jectória com um mínimo de 3 reflexões. 

Para \ip\ > (po o postulado está já demonstrado, porque sabemos, da prova anterior, 
que nesse caso as trajectórias são sempre de 3 reflexões. Faltará então demonstrar para 
—ífo < íf < ipo, mas dada a simetria da cavidade apenas teremos necessidade de provar a 
veracidade da afirmação para o intervalo O < ip < ipo- 

Basearemos a nossa demonstração em algumas das deduções que fizemos na prova anterior 



( ^4.1[ ), sendo-nos especialmente úteis as ilustrações da Figura 12 Assumiremos ainda como 
verdadeira a seguinte premissa: "Se a 2^ reflexão acontecer na mesma face da cavidade onde 
ocorreu a 1^ reflexão, haverá necessariamente uma terceira reflexão". Dispensamo-nos de 
provar este princípio por nos parecer evidente. 

Para O < íp < ípo a primeira reflexão tanto pode ocorrer na face do lado esquerdo como 
na do lado direito. Analisemos cada um dos casos em separado. 

1°' reflexão no lado esquerdo 

Sendo O < tp < ípo poderemos ter as duas primeiras reflexões na face do lado esquerdo, 
estando neste caso garantida, como assumimos atrás, a existência de 3 ou mais reflexões. Caso 
contrário, sendo a segunda reflexão no lado direito, uma parte inicial da trajectória poderá ser 



sempre representada pelas primeiras três ilustrações da Figura 12 (assumindo O < ip < ípo), as 



quais garantem, também nesse caso, a existência de uma terceira reflexão B3. Para provarmos 

> 

o que acabámos de dizer bastar-nos-á demonstrar a natureza ascendente do subtrajecto B1B2. 



Fixemos na parábola do lado esquerdo (ilustração (a) da Figura 12 ) o primeiro ponto de 



reflexão Bi. Para qualquer que seja Bi é sempre possível traçar mos u m subtrajecto inicial 



BqBi com origem num ângulo de entrada (p > ipQ. Como em ^.1 (pág. 16) mostrámos 



ser ascendente o subtrajecto -61-62 que se seguiria, o mesmo necessariamente sucederá para 

> 

qualquer que seja O < (/? < ipQ, dad o qu e neste caso -69-61 apresentará um declive negativo 



mais acentuado. Uma vez que em ^.1 (pág. 19) caracterizámos -62-63 apenas com base na 



natureza ascendente do subtrajecto precedente -61-62, as conclusões a que chegámos para 
.62-63 são igualmente válidas para d < ip < ipQ. 

1°" reflexão no lado direito 



25 



Também neste caso poderemos ter as duas primeiras reflexões na face do lado direito, 
ficando garantida a existência de 3 ou mais reflexões. Se isso não acontecer, teremos nec- 
essariamente uma trajectória com o aspecto da trajectória B0B1B2B3 ilustrada no esquema 



da Figura 15 onde surgem igualmente representados dois trajectos auxiliares (a tracejado), 
^0-^1^2 e A1B2A3, que, ao passarem pelos focos das parábolas, apresentam o subtrajecto 
posterior à reflexão horizontal. Tendo por base as leis de reflexão, podemos sucintamente de- 




Ao L- Ai 




Figura 15: Esquema ilustrativo para estudo da trajectória de partículas com ângulos de 
entrada O < (p < ipo. 



duzir o seguinte: como o ângulo A0B1À2 tem que ser interior ao ângulo BqBiB2, concluímos 
que B1B2 é de natureza ascendente; como B1B2B3 é necessariamente um ângulo interior a 
A1B2A3, concluímos que B^ tem que se situar entre Ai e A^, o que garante a existência da 
terceira reflexão. Fica assim concluída a demonstração do Teorema [2} 



5 Outras possíveis aplicações 

Para além de maximizar a resistência newtoniana, é entusiasmante veriflcar que as potencial- 
idades da forma Dupla Parábola por nós encontrada, podem vir a revelar-se também muito 
interessantes noutros domínios de interesse prático. 

Se revestirmos a parte interior da cavidade Dupla Parábola com uma "superfície" polida, 
a trajectória da luz no seu interior será descrita pelos princípios da óptica geométrica, des- 
ignadamente propagação rectilínea da luz, leis de reflexão e reversibilidade da luz. Assim, 
os modelos computacionais que foram por nós desenvolvidos para simular a dinâmica de bil- 
har no interior da cada uma das formas estudas (onde se consideraram colisões de partículas 
perfeitamente elásticas) são igualmente válidos quando o problema passa a ser de natureza 
óptica. Olhemos estão para a forma 2D por nós encontrada nessa nova perspectiva. 

Dadas as características de reflexão que apresenta a forma Dupla Parábola, rapidamente 
vislumbramos-lhe uma propensão natural para poder vir a ser usada com grande sucesso no 
desenho de retrorreflectores. Retrorreflectores são dispositivos que enviam a luz ou outra 
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radiação incidente de volta à fonte emissora. Idealmente, o retrorreflector deve desempen- 
har essa função independentemente do ângulo de incidência, algo que não acontece com os 
dispositivos existentes. Nesses dispositivos a superfície retrorreflectora é composta por uma 
unidade óptica de reflexão ou, na maior parte dos casos, por um conjunto considerável desses 
elementos em tamanho reduzido. Esses elementos ópticos individuais — responsáveis pela 
inversão do fluxo incidente — podem ser pequenos elementos esféricos com propriedades 
ópticas adequadas, mas o mais habitual é serem pequenos cantos cúbicos retrorreflectores — 
cantos formados por três espelhos mutuamente perpendiculares — , que como se sabe apenas 
invertem a radiação para determinados ângulos de incidência. 

Como vimos, a Dupla Parábola, embora não garanta a inversão perfeita de toda a radiação 
incidente, desempenha essa função com grande sucesso (os pontos (99, (/?"*") concentram-se nas 
proximidades da diagonal 99 = Lp^ , Figura [sj : garante uma grande aproximação das direcções 
dos fluxos incidente e reflectido para uma parte significativa dos ângulos de incidência, e 
mesmo para os restantes não permite que o desfazamento atinja valores elevados. Antevemos, 
por isso, bastante promissora a sua possível utilização na definição de novas geometrias para 
os elementos ópticos que compõem as superfícies retrorreflectoras. 

Os dispositivos retrorreflectores, embora sejam usados nas mais diversificadas áreas tec- 
nológicas, como é o caso por exemplo das comunicações ópticas em espaço aberto, é na 
indústria automóvel e sinalização rodoviária que, sendo utilizados de forma massiva, nos 
apercebemos mais facilmente da sua utilidade. Vejamos, na secção que se segue, uma forma 
possível de tirarmos partido do nosso resultado. 

Retrorreflectores na indústria automóvel e sinalização rodoviária 

Hoje em dia todas as principais vias de circulação automóvel dispõem de dispositivos retror- 
reflectores, complementares à sinalização horizontal, instalados fora da superfície pavimentada 
e a uma altura predeterminada — designados delineadores — , com o objectivo de sinalizar de 
forma clara a geometria da via, mesmo perante condições de visibilidade reduzida como é o 
caso da condução nocturna ou sob condições climatéricas adversas. Cientes do importante pa- 




Figura 16: Exemplo de retrorreflector do tipo delineador (imagem retirada do sítio Web da 
BRISA — Auto-estradas de Portugal). 

pel que esse tipo de sinalização desempenha no aumento dos níveis de segurança rodoviária, 
tentamos, com base nos resultados obtidos, propor um esquema diferente de dispositivos 
retrorreflectores que, reflectindo a luz de forma mais eficiente, apresentarão certamente coe- 
ficientes de intensidade luminosa mais elevados (indicador normalmente usado para medir o 
desempenho dos retrorreflectores). 

Embora a forma por nós encontrada tenha existência no espaço bidimensional, é possível, 
para determinadas aplicações práticas, construir a partir dela formas 3D que mantenham 
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as propriedades ópticas adequadas. Vejamos o caso dos retrorreflectores usados como delin- 
eadores nas vias de comunicação. Como a altura a que se encontram do pavimento, tanto 
os olhos do condutor (ponto de observação) como os faróis do veículo (fonte de luz), é sen- 
sivelmente a mesma para a generalidade dos automóveis, pode-se facilmente projectar um 
reflector que garanta que o raio de luz reflectido, quando projectado no plano vertical que 
une o veículo ao reflector, passe ao nível do ponto de observação. Esse efeito é conseguido 
com espelhos de superfície vertical colocados a uma altura equivalente à média das alturas a 
que se encontram o ponto de observação e a fonte de luz — aproximadamente 75 cm acima 
do pavimento da faixa de rodagem — , tal como ilustrado na figura [17} Assim, a parte da 




Figura 17: Esquema ilustrativo da luz reflectida por dispositivos com superfície reflectora 



geometria da superfície reflectora que falta ainda definir é neste caso a forma da sua projecção 
no plano horizontal. Será essa forma que determinará o grau de aproximação dos raios de 
luz incidente e reflectido, quando projectados num plano horizontal, e que idealmente de- 
verão ter a mesma direcção — ver figura IT8| Tratando-se de uma curva 2D, a forma Dupla 



"""""""""■^^^^^^^^^^^^ retrorreflector 

berma 

Figura 18: Esquema ilustrativo da luz reflectida por retrorreflectores ideais. 



Parábola será seguramente a melhor escolha, pois é aquela que em nosso entender apresenta 
melhor desempenho — ver figura 19 Assim, a superfície retrorreflectora 3D que propomos é 



Retrorreflector 
proposto 



berma 




Figura 19: Esquema ilustrativo da luz reflectida pelo retrorreflector proposto. 



a parte interior da superfície varrida pelo deslocamento vertical do contorno da nossa forma 
2D contida num plano horizontal, tal como se ilustra na figura 20 Um inovador dispositivo 
retrorreflector pode então ser construído integrando o elemento proposto como unidade única 
de reflexão ou um conjunto desses elementos distribuídos numa área rectangular segundo uma 
disposição matricial. 

Se se entender que também no plano vertical o dispositivo deve reflectir a luz de forma a ter 
a mesma direcção da luz incidente, a solução passará por se introduzir uma pequena alteração 
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Figura 20: Forma 3D da superfície proposta para elemento retrorreflector. 



na forma do elemento retrorreflector: para além da superfície vertical, a forma passará a ser 
delimitada no topo e na sua base por duas superfícies planas horizontais. Sabemos que, 
com esta configuração, o sentido da luz é invertido sempre que ocorra um número ímpar de 
reflexões nas superfícies horizontais — lados superior e inferior. Isto garante-nos que uma 
boa parte da luz reflectida toma efectivamente a direcção da luz incidente. Em [7] foram já 
realizadas algumas simulações numéricas para o caso 3D, tendo os resultados apontado para 
o bom desempenho desta configuração — trata-se da forma 3D de maior resistência que se 
conseguiu obter. 

O mesmo tipo de solução pode ser usado para os retrorreflectores instalados na parte 
posterior das viaturas. Também neste caso o reflector encontra-se normalmente a uma altura 



intermédia entre a fonte de luz e o ponto de observação, como se ilustra na figura 21 



O: 



Figura 21: Esquema da reflexão de luz entre dois veículos. 



Poderão também ser conseguidos óptimos resultados com o uso combinado do nosso ele- 
mento reflector com outros reflectores já existentes, dado que se podem complementar. Por 
exemplo, as formas em ângulo recto têm um óptimo desempenho para ângulos de incidência 
baixos, mas deixam de ter qualquer eficácia para ângulos superiores a 45°; pelo contrário, o 
nosso elemento reflector tem um desempenho irrepreensível para ângulos superiores a 45°, e 
perde alguma eficácia em ângulos de incidência de baixa amplitude. 

Ainda que as reais mais-valias do nosso modelo de retrorreflexão apenas possam ser ver- 
dadeiramente comprovadas através da realização de ensaios em laboratórios especializados 
que permitam a análise comparativa com outros modelos, estamos em todo o caso seguros da 
sua eficácia como retrorreflector, em situações em que os emissores/receptores se movimentem 
num mesmo plano, ou então, se quisermos ser mais precisos, em todas as situações em que 
seja possível encontrar um plano bidimensional que intercepte perpendicularmente e a meio 
o eixo que une o emissor ao receptor de todos os pares emissor /receptor admissíveis. Nos 
exemplos que demos considerámos ser essa a situação, embora saibamos que na prática não 
se passe exactamente assim. 
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6 Conclusão 



Na continuação do estudo realizado anteriormente pelos autores em |181ll9j. conseguiu-se, com 
o trabalho agora apresentado, obter um resultado original que nos parece de grande alcance: os 
algoritmos de optimização convergiram para uma forma geométrica muito próxima da forma 
ideal — a Dupla Parábola. Trata-se de uma forma de rugosidade que confere uma resistência 
quase máxima (muito próxima do majorante teórico) a um disco que, para além de se deslocar 



num movimento translacional, rode lentamente sobre si mesmo. Na Figura 22 é mostrado um 
desses corpos. Atendendo a que o contorno do corpo apresentado é integralmente formado 




Figura 22: Corpo 2D (quase) óptimo. 



por 42 cavidades O, com a forma da Dupla Parábola, cada uma das quais com uma resistência 
relativa de 1.49650, de (|2| e (|4| concluímos que R(B) = ''"]j/j^^ .R(0) « 1.4951 é a resistência 
total do corpo, um valor 49.51% acima do valor de resistência do disco de contorno liso 
correspondente (o menor disco que inclua o corpo). Sabemos que se o corpo for formado por 
um número suficientemente elevado de pequenas dessas cavidades, a sua resistência atingirá 
mesmo o valor 1.4965, mas o exemplo apresentado é suficiente para percebermos o quão 
próximos ficámos do conhecido majorante teórico (50%). 

Ainda que o valor da resistência da Dupla Parábola tenha sido determinado numerica- 
mente, o estudo que se realizou em Q com o objectivo de se caracterizar analiticamente 
as reflexões no interior dessa forma de cavidade, consolida os resultados apresentados. Con- 
seguimos demonstrar algumas propriedades importantes que ajudam a perceber o elevado 
valor de resistência que se obteve. Tentaremos no futuro desenvolver outros estudos teóricos 
que venham a permitir consolidar ainda mais este resultado. Por exemplo, um problema 
interessante em aberto consiste em delimitar o desfazamento entre os ângulos de entrada e 
de saída para as trajectórias de 3 reflexões — para as restantes (trajectórias com 4 ou mais 
reflexões) sabemos já que \ip — ip~^\ < 2lpq ~ 2 x 19.47°. 

A Dupla Parábola é efectivamente um resultado de grande alcance prático. Para além de 
permitir elevar a resistência newtoniana quase ao seu majorante teórico, essa forma geométrica 
revela-se também muito interessante noutros domínios de aplicação. Dadas as suas carac- 
terísticas de reflexão, rapidamente vislumbramos-lhe uma propensão natural para poder vir 
a ser usada com grande sucesso no desenho de retrorreflectores — veja- se em ^ o estudo 
exploratório sobre a sua possível utilização na sinalização rodoviária e indústria automóvel. 
Mas a sua utilidade estende-se a muitos outros domínios de aplicação importantes, como será 
certamente o caso dos sistemas de comunicações ópticas em espaço aberto ou dos sistemas 
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ópticos que integram alguns dos equipamentos tecnológicos ligados às ciências da saúde e à 
indústria em geral. Na continuação deste trabalho, tencionamos vir também a desenvolver 
investigações que permitam a análise comparativa de superfícies baseadas na Dupla Parábola 
com outros modelos de retrorreflexão existentes. Nesse estudo usaremos, em substituição da 
função resistência, uma função de custo mais adequada à avaliação das capacidades retror- 
reflectivas — uma função que contabilize devidamente o desfazamento angular entre os fluxos 
incidente e reflectido. 

Uma incursão no domínio tridimensional, realizada em ^1, mostrou também ela, ser a 
Dupla Parábola uma forma de cavidade muito especial. A nossa convicção no seu virtuosismo 
saiu realmente reforçada quando obtivemos o melhor resultado para o caso 3D. Esse resultado 
foi conseguido com uma cavidade cuja superfície é a área varrida pelo deslocamento da curva 
Dupla Parábola na direcção perpendicular ao seu plano. Tendo o valor da sua resistência 
(R = 1.80) ficado um pouco aquém do majorante teórico para o caso 3D {R = 2), superar 
esse valor será também um desafio interessante a considerar futuramente. 

Já para o caso 2D pressentimos maior dificuldade em virmos a superar o resultado já 
alcançado — quer pela proximidade a que este se encontra do majorante teórico, quer pelo 
facto de termos já realizado, sem sucesso, uma série de investigações visando esse objectivo. 
Tendo ainda em conta que o majorante 1.5 apenas significa a não existência de formas que 
superem esse valor de resistência, a forma por nós encontrada poderá mesmo tratar-se de 
uma solução óptima. A confirmar-se esta hipótese — o valor da resistência que caracteriza 
a Dupla Parábola tratar-se efectivamente do limite máximo que é possível alcançar com 
formas bidimensionais reais — o nosso resultado ganhará ainda uma importância redobrada. 
Embora não se perspectivem fáceis desenvolvimentos, este é um importante problema que fica 
em aberto à espera de futuras contribuições que, se não vierem a superar o nosso resultado, 
venham a permitir reforçar a nossa conjectura. 
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